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ABSTRACT 


For  a  stochastic  epidemic  of  the  type  considered  by  Bailey  [1]  and 
Kendall  [3],  Daniels  [2]  showed  that  "when  the  threshold  is  large  but 
the  population  size  is  much  larger,  the  distribution  of  the  number  remain¬ 
ing  uninfected  in  a  large  epidemic  has  approximately  the  Poisson  form." 

A  simple,  intuitive  proof  is  given  for  this  result  without  use  of  Daniels' 
assumption  that  the  original  number  of  infectives  is  "small".  The  proof 
is  based  on  a  construction  of  the  epidemic  process  which  is  more  explicit 
than  the  usual  description. 


Key  words:  limiting  Poisson  distribution,  stochastic  epidemic 


eor* 


V  ,.T  « CAikI 
!  j  TAB 

.  .c.'on«ad 


'  lr;' ntl on/  — 

«  r*f  £  ,  5P 


.  .r  -n  al 


ii 


ON  THE  ASYMPTOTIC  DISTRIBUTION  OF  THE  SIZE  OF 


A  STOCHASTIC  EPIDEMIC 

1.  Introduction 

Consider  a  population  which  at  time  t  ■  0  consists  of  X(0)  -  n 
healthy  individuals  and  Y(0)  “  ■  individuals  with  a  contagious  infection. 

An  epideaic  in  such  a  population  is  often  modeled  by  a  continuous- time  Markov 
process  as  follows.  (See,  for  example,  Bailey  [1]  or  Kendall  [3].)  If 
X(t)  and  Y(t)  are  the  numbers  of  healthy  individuals  and  infectious 
individuals,  respectively,  present  at  tine  t,  then  the  transition 
probabilities  are  given  by 

P{(X(t+6),Y(t+6))  -  (x',y*)|(X(t),Y(t))  -  <x,y)} 

-  xyd  +  o(6)  for  (x’,y')  -  (x-l,yfl) 

-  pyfi  +  o(6)  for  (x',y')  -  (x,y-l) 

-  l-xy6-py6  +  o(6)  for  (x'.y’)  -  (x,y)  . 

The  transitions  listed  represent  the  infection  of  a  healthy  individual, 

the  removal  of  an  infectious  individual  from  the  population,  and  "no  change", 

respectively.  All  other  possible  transitions  in  [t,t+6]  are  assumed  to 

have  collective  probability  o(6) .  The  positive  constant  p  is  variously 

called  the  threshold  or  the  relative  removal  rate  of  the  epidemic. 

Note  that  states  of  the  form  (x,0)  are  absorbing,  so  that  no  more 

transitions  occur  after  the  last  infectious  individual  has  bean  removed. 

Absorption  at  (x,0)  means  that  x  individuals  have  escaped  infection 

at  the  end  of  the  epidemic  and  that  X(»)  •  x,  Where  X(«)  »  11m  X(t) .  It 

t  ♦  • 

is  easy  to  see  that,  with  probability  one,  acme  such  absorbing  state  la 


reached  eventually.  Most  of  the  work  on  this  stochastic  epidemic  aodel 
has  been  directed  toward  finding  the  distribution  of  the  number  of  new 
infections  occurring  during  the  course  of  the  epidemic.  This  is  equivalent 
to  finding  the  distribution  of  the  absorbing  state.  See  Nagaev  and  Startsev 
[4]  for  an  asyaptotlc  analysis  of  Halting  cases  other  than  the  one 
considered  in  this  paper. 

Under  the  assumption  that  Y(0)  ■  a  is  "small",  Daniels  [2]  shows 
that  "when  the  threshold  is  large  but  the  population  size  is  much  larger, 
the  distribution  of  the  number  remaining  uninfected  in  a  large  epidemic 
has  approximately  the  Poisson  form  with  deterministic  mean  ne~n^."  However 
the  asymptotic  approximations  used  by  Daniels  in  obtaining  this  result  are 
not  very  enlightening t  and  he  speculates  that  "there  must  be  a  direct  argu¬ 
ment  in  terms  of  the  epidemic  process  itself  to  explain  this  Polssonllke 
behavior . " 

This  paper  gives  a  rather  intuitive  proof  of  Daniels'  result.  In 
addition,  the  assumption  that  m  is  small  is  dispensed  with.  The  proof  is 
based  on  a  construction  of  the  epidemic  process  which  is  more  explicit  than 
the  usual  description  given  above. 

2.  Construction  of  the  Epidemic  Process 

Let  the  n  originally  healthy  individuals  be  Indexed  by  i,  1  £  i  £  n, 
and  let  the  m  originally  infectious  individuals  be  indexed  by  J, 

1  £  J  £  m.  Let  and  be  i.i.d.  random  variables  with 

density  pe“pt  on  [0,®) .  Individual  j  in  the  original  infectious  group 
will  remain  infectious  for  time  units  before  removal  from  the  population 

Individual  i  in  the  original  healthy  group  will  remain  infectious  for  r^ 
time  units  if  he  becomes  infected. 
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Let  b«  1.1, d.  random  variables  with  density  e-t  on 

(0,«»)»  Independent  of  the  r^’a  sad  r^'s.  The  vsrlsble  will  be 
thought  of  ss  the  "resistance  to  Infection"  of  Individual  1  In  the 
original  healthy  group.  Let  be  the  associated  order  statistics, 

so  that  <  ...  <  Let  r'k'  *  rA  if 

Now  let  the  epidealc  process  proceed  as  follows.  The  originally 
Infected  Individual  j  remains  in  the  population  for  tine  units, 

after  which  he  is  removed.  The  healthy  individual  1  accumulates  "exposure 
to  infection”  at  a  rate  equal  to  the  number  of  infected  individuals  present. 
When  the  total  exposure  to  infection  of  healthy  individual  1  reaches  1^, 
individual  i  becomes  infected  and  then  remains  in  the  population  for  an 
additional  time  units  before  removal. 

It  remains  to  be  shown  that  the  resulting  process  is  Markov  with 
the  correct  transition  probabilities.  Suppose  (X(t),Y(t))  *  (x,y).  The 
probability  that  a  particular  infected  individual  is  removed  in  the  time 

A 

interval  [t,t+6]  is  p6  +  o(6)  because  the  distribution  of  the  r^'s  and 
r^'s  has  constant  hazard  rate  p.  The  probability  that  exactly  one  of 
the  y  infected  individuals  is  removed  in  [t,t+6]  is  therefore  pyd  +  o(5). 
The  probability  that  a  particular  one  of  the  x  healthy  individuals  will 
become  infected  in  [t,t+6]  is  yd  +  o(6),  so  that  the  probability  that 
exactly  one  of  the  healthy  individuals  becomes  infected  is  xyd  +  o(6).  It 
follows  that  the  transition  probabilities  are  as  desired.  The  Markov  property 
follows  from  the  memory less  property  of  the  exponential  distribution. 

Let  v  be  the  number  of  new  infections  occurring  during  the  course  of 

m  A 

the  epidemic.  If  I...  >  Jr.,  then  all  originally  infectious  individuals 

W  j»l  j 

are  removed  before  the  resistance  to  infection  of  any  healthy  individual  has 
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been  exceeded,  so  that  v  -  0.  Otherwise,  the  originally  healthy 
individual  associated  with  becomes  the  first  new  Infection,  and 

v  >  1.  An  easy  Induction  argument  shows  that  V  +  1  is  the  smallest  k, 
1  <  k  £  n,  for  which 


m  A  k-1 

A(k)  >  ^  r1  +  ^  r 
w  j-1  J  i-1 


(1) 


If  this  Inequality  does  not  hold  for  any  k,  1  £  k  £  n,  then  v  ■  n. 
Define 

R-  i;,+  l  r(1>  . 

J-1  5  1-1 


Then  R  Is  the  amount  of  "exposure  to  infection"  withstood  by  those 
individuals  who  remain  healthy  at  the  end  of  the  epidemic,  and  X(«)  ■  n  -  v 
is  the  number  of  1^'s  greater  than  R. 

3.  Statement  of  the  Theorem  and  an  Outline  of  the  Proof 

Consider  a  sequence  of  epidemics  with  parameters  n^,  m^,  and  p^, 

1  £  k  <  «».  Let  be  the  number  of  new  infections  in  a  realisation  of 
the  k-th  epidemic,  so  that  X^(®)  -  is  the  number  of  individuals 

who  escape  infection. 


Theorem.  If  %  ■*“  **  P*  ■*’  *»  and 


}  b  ,  0  <  b  <  •  , 


then  Xj^W  converges  in  distribution  to  a  Poisson  random  variable  with 


b. 
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Suppressing  the  subscript  k,  we  can  summarize  the  proof  as  follows. 

The  fact  that  p  is  o(n+m)  Implies  that,  with  high  probability,  all  but 

a  tiny  fraction  of  the  population  becomes  infected.  Thus,  R  will  be  close 
m  n 

to  J  r.  +  Jr.,  which  is  in  turn  close  to  (m+n)/p  with  high  probability. 
j-1  2  i-1 

If  R  is  close  enough  to  <m+n)/p,  then  X(®)  will  equal  the  number  of  J^'s 
which  are  greater  than  (rrhn)/p.  The  number  of  I^'s  greather  than  (mfn)/p  has 

the  distribution  Binomial  (n,e  t  which  converges  in  distribution  to 

Poisson  (b). 

4.  Proof  of  the  Theorem 

The  subscript  k  will  again  be  suppressed. 

Taking  logarithms  in  ne  ^m+n^P  -+■  b  yields 

log  n  -  ♦  log  b  , 

so  that  p  *V»  Thus,  p  is  o(m+n),  but  (mfn)^  is  o(p)  for 

0  <  y  <  1. 

2 

Lemma  1.  Let  0  <  e  <  1  be  given.  If  p<  e  (m+n),  then 
P{X(<»)  >  e(m+n)}  <  <  em  <  e  . 

Proof  of  Lema  1. 

.4s  in  Kendall  [3],  we  can  view  the  population  of  infected  individuals 
as  a  continuous-time  birth-and-death  process  with  a  variable  birth-rate. 

The  ratio  of  death-rate  to  birth-rate  is  p/X(t),  which  is  less  than 
p/(e(*fn)) ,  until  X(t)  <  e(m+n).  The  probability  that  a  birth-and-death 
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process  starting  at  m  and  with  a  death-rate  to  birth-rate  ratio  q  <  1 
is  ever  absorbed  at  zero  is  qQ. 

Leona  2,  Let  0  <  e  <  1.  For  n  sufficiently  large* 

P{R  <  .<L-2eK»*g>}  <  2e  . 

Proof  of  Lemma  2. 

By  Lenma  1,  R  is  greater  than  the  sum  of  the  first  (1-e) (m+n) 

A  A  /v  /IN  O) 

terns  of  iri»r2* . . . ,rm,  r'  ,rv  ,...,r  }  with  probability  greater 

than  (1-e)  for  sufficiently  large  n.  The  sum  of  the  first  (1-e) (m+n) 
terms  has  mean  (1-e) (mfn)/p  and  variance  (1-e) (m+n) /p  .  The  Chebyshev 
inequality  now  Implies  the  lemma. 

Lemma  3.  Let  0  <  e  <  1.  For  n  sufficiently  large, 

P{X(»)  >  n3€)  <  3e  . 

Proof  of  Lemma  3. 

The  number  of  Jl^'s  which  are  greater  than  (l-2e)  (m+n)/p  is 
Binomial  (n,e  ^  ^e)  (®+n) / .  This  distribution  has  mean 

M-(l-2€K»h,)/p  n(t)a-2e)  .  bl-2£n2c  _ 


Application  of  Lenma  2  and  the  Markov  inequality  finishes  the  proof  of 

Prom  Leema  3,  it  is  easy  to  see  that  X(«)  is  o(./o)  in  probability 
Thus,  except  on  a  set  of  small  probability,  R  is  greater  than  the  sun  of 
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